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We examine the spin and charge excitations in antiferromagnetic iron pnictides by mean-field 
calculations with a random phase approximation in a five-band itinerant model. The calculated 
excitation spectra reproduce well spin-wave dispersions observed in inelastic neutron scattering, 
with a realistic magnetic moment for CaFe2As2. A particle-hole gap is found to be crucial to 
obtain consistent results; we predict the spin wave in LaFeAsO disappears at a lower energy than 
in CaFe2As2. We analyze that the charge dynamics to make predictions for resonant inelastic x-ray 
scattering spectra. 
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I. INTRODUCTION 

Iron pnictide superconductors have been intensely 
studied since the discovery of superconductivity in 
LaFeAsOi-zFa—, and achieving critical temperatures 
(T c ) over 50 K. Such a high T c evokes the high-T c 
cuprates, and the presence of the antiferromagnetic 
(AFM) order in the parent compounds suggests a close 
connection between the superconductivity and the mag- 
netic order in iron pnictides as well. From this similarity, 
one may expect to elucidate the origin of the high-T c 
mechanism, and the relation with the magnetic proper- 
ties of these materials, in the same framework. 

However, there exists one crucial difference between 
these systems — iron pnictides are not a strongly corre- 
lated system such as the cuprates. Although the parent 
compounds of the high T c cuprates are a Mott insula- 
tor, those of the iron pnictides are a metal. Because 
of the itinerant feature in iron pnictides, the magnetic 
properties are not so readily understood as those in the 
insulating AFM systems. 

The most powerful tool to observe the spin excitations 
is inelastic neutron scattering, and the observations in 
CaFe 2 As 2 and BaFe2As2 (the so-called 122 system) have 
discovered the following characteristics of the spin exci- 
tation. (I) The spin wave excitation is extended up to 
a high energy (not less than 100 meV) 2 -^. (II) The ob- 
served spin wave excitation is anisotropic in the planed; 
this anisotropy has been also found in the paramagnetic 
(PM) state in the electron-doped system-, (III) The 
wave-vector-averaged intensity of the excitation spectra 
is lower in the AFM case than that in the PM case below 
80 meV— . Most of the experimental data have been ana- 
lyzed in terms of a local-spin model^, except for RefA 

Some other properties of the spin excitations in AFM 
metals have been illuminated in theoretical work££, but 
a sufficient understanding of the itinerant AFM phase is 
not yet provided. Since understanding the itinerant AFM 
system is a key to progress in the study of iron pnictide 
superconductors, we aim to find an appropriate recipe 
that can describe the basic properties of the AFM phase. 
Our recent work 9 on modeling the AFM phase of the 



iron pnictides revealed the weakly ordered character of 
the ground state by analysis of the optical conductivity. 
Based on this success of the ground-state analysis, we 
attempt here to analyze the spin excitations in the AFM 
ground state by extending the mean-field approach. 

In the next section, we will show our mean-field cal- 
culation and the results obtained from those calcula- 
tions. In the section after the next, we investigate the 
spin excitations in an itinerant AFM system by mean- 
field calculations with a random phase approximation 
(RPA) within the five-band model that can well de- 
scribe iron pnictides. By comparing our results with 
those of the experiments for the 122 system, we demon- 
strate how the mean-field model consistently describes 
the spin-excitation properties, in terms of the presence of 
the spin wave excitation up to a high energy [character- 
istics (I)], its anisotropic behavior [characteristics (II)], 
and the spectral intensity relation between the PM and 
AFM states [characteristics (III)] . Based on the calcula- 
tion of the particle-hole excitation spectra, we show that 
such excitations cause damping of the spin wave excita- 
tion, predicting that the spin wave excitation in LaFeAsO 
(the so-called 1111 system) disappears in a lower energy 
than in the 122 case. In the following section, we discuss 
the difference between the transverse and longitudinal 
modes of the spin excitations, and also provide a pre- 
diction about a charge excitation for a resonant inelastic 
x-ray scattering (PJXS) experiment. 

II. MEAN-FIELD FIVE-BAND MODEL 

Considering an Fe square lattice (the Fe-Fe bond 
length ao set to be unity; the x and y directions, along 
the nearest Fe-Fe bonds), we start with the five-band 
mean-field Hamiltonian represented by the ordering vec- 
tor Q = (|rj,0) (Nq — 2 for AFM, and 1 for PM): 
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where c k v creates an electron with wave vector k and 
spin a at orbital fi. The diagonal component of Hmf IS 

H^(k,a) = J2^x,A y ;^iy)e^ h+lQ > A + (2) 

A 

where t(A x , A y ; /i, v) and e M are the tight-binding ener- 
gies presented in RefJ^, and A = (A x ,A y ). The off- 
diagonal component is 
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where U is the intraorbital Coulomb interaction, J is the 
Hund coupling, and the pair hopping is set equal to J. 
Hamiltonian ([T]) is derived from the tight-binding+([7, J) 
Hamiltonian represented in Rcf.^ by retaining the spin- 
density- wave order parameter defined as 
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where N is the number of k points in the first Brillouin 
zone (BZ) of the five-band PM system, and I ^ 0. 

To obtain the ground state, we solve mean-field equa- 
tions self-consistently to obtain the quasiparticle state 

llns =^2'>P^l;n( k ^) C k+lQ fl a ( 5 ) 
//./ 

with the energy Ek,n,a- The ground state is represented 
as a set of quasiparticles distributed according to the 
Fermi distribution function /. Different parameter sets 
yield different strength of the order. We evaluate the 
order strength from the magnetic moment: 
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The parameter set U = 1.1 eV and J — 0.2 eV yields 
M = 0.4/ib corresponding to the 1111 system^ (we refer 
to this as the 1111 model). The Fermi surface of this 
model is plotted in Ref£. To simulate the 122 system 1 ^, 
we use the 10% larger U = 1.2 and J = 0.22 13 , which 
yields M = 0.8/ie (the 122 model). The symmetry- 
broken Fermi surface of this model is plotted in Fig. Q] 
together with the symmetric one in the PM case. We 
note that Dirac cones appear near (0, 0) along the k x 
axis as a small electron pocket in the AFM case and 
that the presence of the Dirac cones affects the trans- 
port properties 14 . The order strength is also evaluated 
from the partially opened gap, which is estimated to be 
- 0.2 eV (122) and - 0.1 eV (1111) from the density of 
states (for the 1111 model, shown in Rcf.-). The gap in 
the single-particle excitation near the T point opens at 
around k = (±0.27r, ±0.27r), as can be found in Fig. [T] 
From the orbital-resolved density of states plotted in 
Fig. [3J it is also found that the xy orbital component 
is dominant near above the Fermi level. 





FIG. 1. (Color online) Fermi surfaces for the PM case (a) and 
the 122 case (b). 



III. DYNAMICAL SUSCEPTIBILITIES 

To investigate the spin excitations around the mean- 
field ground state, we calculate the dynamical susceptibil- 
ity in the spin transverse channel by RPA in Matsubara 
form, 

xtx (fei , k 2 , iw) = xt (fei , k 2 , iw) 
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where xo is the bare susceptibility whose explicit form 
is given below, the superscript + represents a pair of a 
down-spin Q.) hole and an up-spin (t) electron, and the 
nonzero elements of the interaction matrix V ++ arc 
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The dynamical susceptibilities of the longitudinal 
modes are calculated from 
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where the orbital indices are omitted together with their 
summations, which are taken in the same manner as in 
the transverse case, Eq. ([7]). The nonzero elements of the 
interaction matrix V" are 



U for A = r = fi 

J for A = t^/i 
U - 3 J - J8 a . a 

,J-(U- 2J)8 a ^ for A = v ^ fi = t 



for A = fi T = v v 7 



The bare susceptibility in Matsubara form is repre- 
sented with the wave functions and the quasiparticle en- 
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FIG. 2. (Color online) The orbital-resolved density of states 
of the majority spin (a) and the minority spin (b). 
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where the set of spins (a, a') takes (f, t)i (4-i l)j an d (t\ 4) 
for s =f, 4-, and +, respectively. The summation of p 
runs over the reduced magnetic BZ. 

We evaluate the imaginary part of the dynamical and 
the bare susceptibilities 



X 
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to determine the collective and the individual excitations, 
respectively, where rj is set to 0.01 eV. 
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FIG. 3. (Color online) ( X + )" [(a)-(c)] and [(d)-(f)] for 

the M = 0.8hb [(a), (d)], the M = 0.4/i B [(b), (e)], and the 
M — O.OfiB [(c), (f)] cases. The color bars are common in (a)- 
(c) and in (d)-(f), respectively. Intensities higher than 200 are 
not taken account of in (a) and (b) for better visibility: The 
maximum intensities are ~ 2800 (a) and ~ 1000 (b). 



IV. SPIN TRANSVERSE EXCITATIONS 

Figures Ela) and 02[b) show the transverse spin- wave 
excitation spectra in the 122 and 1111 models, respec- 
tively. We find that a collective mode appears at (n, 0) 
and persists up to ~ 0.2 eV in the 122 model, reproducing 
the characteristics (I). Similar spin- wave excitations were 
obtained within an effective three-band model 15 . Above 
~ 0.2 eV, this spin wave excitation is damped; this damp- 
ing feature has also been noted in Refj£. The 1111 model 
shows the damping at lower energy (~ 0.1 eV). 

A possible cause of this damping is particle-hole exci- 
tations, as mentioned in Ref. 2 ; however, there remains 
a question of why only the high-energy excitations are 
damped despite the itinerant system. To investigate this, 
we plot the particle-hole excitation spectra in Figs. El[d) 
andfSJe). The (Xo)" spectra exhibit the strong excitation 
spectra above 0.2 eV (0.1 eV) in the 122 (1111) model. 
This characteristic energy is consistent with the damping 
feature in (x + )"; therefore, it is clear that the damping 
is caused by the particle-hole excitations. Notice that 
the threshold energies, above which the particle-hole ex- 
citations occur, correspond to the partially opened gap 
estimated from the density of states. 

For the PM case, there appears a broad excitation 
structure around (it, 0) in (x + )" [Fig. El (c)]. This broad 
structure is caused by the particle-hole excitations — the 
(Xo )" s P ec tra exhibit the gapless feature with some in- 
tensities in the low-energy region [Fig. [3] (f )] . 

The energy dependence of the spin-wave cone is plot- 
ted in Figs. IHe)-U][h). Our results reproduce the neu- 
tron experiments, including the anisotropic structure of 
the spin- wave cone (II). In the PM case [Figs. BJa)-|3|d)], 
on the other hand, the excitations appear in a wide re- 
gion around (tt, 0) at each energy. At low energy [Fig. [U 
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FIG. 5. (Color online) (a) (x™)" for the 122 case, (b) ( X z )" 
for the 122 case, and (c) (x™)" for the PM case, (d) The 
difference between the AFM and PM cases [(c) subtracted 
from (a)], (e) Detail of the squared region in (d). 
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FIG. 4. (Color online) The intensity map of (x + )"(k,uj) for 
the PM [(a)-(d)] and the AFM (M = 0.8fi B ) [(e)-(h)] cases at 
the energy u = 200 meV [(a), (e)], 100 meV [(b), (f)], 50 meV 
[(c), (g)], and 10 meV [(d), (h)]. Intensities higer than 200 
are not taken account of in (g) and (h) for better visibility, 
and the color bar is common in (a)-(h). Plotted in (i) are the 
integrated intensities 4 ^2 k (x + )" f° r the PM (circle) and the 
AFM (square) cases. 



AFM states. The problem with the ringlike structure in 
the PM case requires a further study including such ef- 
fects to be addressed, and this should be discussed in 
future work. Nevertheless, the main point of our in- 
terest, the anisotropic feature, is well reproduced. The 
anisotropy in the excitation spectra of the AFM and PM 
states is owing to the anisotropic structure of the bare 
susceptibility arising from the energy band structure. 

Comparing integrated intensity over the momentum 
space in the AFM and PM cases [Fig. H(i)], we find that 
the PM state admits more spin transverse excitations 
than the AFM state. This result matches the spectral 
intensity relation (III) between the PM and AFM states. 

In our results of the AFM case, no gapped feature 
is found at ~ 7 meV, inconsistent with the experi- 
ments^^— . To reproduce the spin-gap feature, we pre- 
sumably need to take into account effects not included in 
our calculation, such as the single-ion anisotropy 3,6 . 



(d)], strong intensities lie along the ring around (ir, 0) re- 
flecting the imperfect nesting of the Fermi surface — this 
excitation ring is anisotropic. Experimentally observed 
spectra in the PM stated, however, show not such a ring 
structure but a broad spot around (7r,0). This broad 
structure may come from finite-temperature effects that 
are not included in the calculation. Strictly speaking, we 
also need to consider a state-dependent r\ taking into ac- 
count the different scattering properties in the PM and 



V. SPIN LONGITUDINAL AND CHARGE 
EXCITATIONS 

Now we analyze the spin longitudinal and charge ex- 
citation modes. The excitation spectra for the modes 
[x n — X^ + X^i X z = \{X^ — X^)} m the 122 case 
are plotted in Figs. [Sfa) and E^b). The intensity of 
these modes, in contrast to the spin transverse mode, 
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is weak in the low-energy region; this means that the 
particle-hole excitations without a spin flip inside the gap 
is strong enough to damp the low-energy excitations of 
these modes. 

To investigate the charge dynamics, we propose a RIXS 
experiment on these materials, which can measure the 
momentum-resolved charge excitations. For this pur- 
pose, we examine the details of the charge excitation 
mode from a theoretical point of view. Compared to the 
PM state [Fig. [5] (c)], which shows a rodlike structure 
rising from (0,0), the excitations in the AFM state are 
rather broad. This difference arises from the presence 
of the magnetic order with which the charge collective 
excitations cost more energy accompanied by the spin 
excitations of the longitudinal mode. The difference in 
the spectra becomes clearer in the plot of the subtraction 
[Fig. [5] (d)] , where the rodlike excitation structure in the 
PM state becomes narrow and enhanced; in RIXS, this 
change should be observable — for example, the change 
comparable to 40% of the peak intensity is found at 0.4 
eV. 

We also find the difference in the low-energy excita- 
tions around (0,0) [Fig. [S] (e)]. Since the spectra of 
the AFM case shows a broad structure there, low-energy 
charge fluctuations can occur with various wave vectors 
away from (0,0). These low-energy charge fluctuations 
may be related to nematic charge structures observed 
by spectroscopic imaging-scanning tunnel microscopy 19 , 
where an eight-site periodic structure represented by the 
wave vector (ir/A, 0) is observed. 

VI. CONCLUTIONS 

In summary, we have investigated the spin excitations 
in an itinerant AFM system by mean-field calculations 
within an RPA in a five-band model. Our results have 
reproduced the characteristics (I)- (III) observed in ex- 
periments for the 122 system. In the spin-wave excita- 
tion, the states outside the partially opened gap plays an 
important role. Particle-hole excitations across the gap 
with a spin flip cause the damping of the spin- wave exci- 
tations above the threshold energy corresponding to the 
partially opened gap. On the other hand, such excita- 
tions within the gap are too weak to damp the collective 



excitation. We predict that the spin-wave excitation in 
the 1111 system disappears at a lower energy than in the 
122 case since the magnetic moment — directly related 
to the partially opened gap — is smaller in the 1111 sys- 
tem. So far, the inelastic neutron scattering in the 1111 
system has been performed only on a powder sample^; 
the experiment on a single crystal is desired. 

In contrast to the spin transverse mode, the charge and 
spin longitudinal modes are weak even in the low energy 
region, because the particle-hole excitations arising from 
the states inside the gap damp these excitations. From 
the analysis of the charge excitation, we provide a predic- 
tion for a RIXS experiment. A rodlike structure should 
be observed in the difference of the excitation spectra for 
the AFM and PM cases. In addition, we have found that 
the magnetic order causes the charge fluctuations. These 
charge fluctuations lying away from (0, 0) may be related 
to the nematic charge structure observed recently, and 
this structure would not exist in the PM state, where 
such charge fluctuations are found to be weak. This 
structure formation should involve the xy orbital com- 
ponent, which is dominant in the states in the energy 
range 0-50 meV. 

In conclusion, all the above features of the spin and 
charge dynamics are ruled by the magnetic order, and 
the characteristic energies for damping and the strength 
of the charge fluctuations are expected to scale with the 
magnetic order strength, i.e., the magnetic moment. 
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